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We calculate the shot noise power and the Fano factor of a spinless resonant level model. The 
Coulomb interaction which in this model acts between the lead electron and the impurity is consid- 
ered in the first order approximation. The logarithmic divergencies which appeared in the expres- 
sions for shot noise and the transport current are removed by renormalization group analysis. It is 
important, that unlike a few recent statements, Keldysh technique gives an adequate description of 
perturbation theory results. By passing to the bosonized form of the resonance model we show that 
in the strong interaction limit the tunnelling becomes irrelevant and decreases. 

PACS numbers: 72.10.Fk, 72.15.Qm, 73.63.Kv 



Introduction: During the last two decades the trans- 
port properties of quantum impurity have attracted great 
interest as a basic problems in nanophysics [l[ . Recently, 
spinless interacting resonance level model (IRLM) [2| has 
become a subject of a special attention [H, 0, 0,13, Sj- 
This interest, in particular, was related to the extend- 
ing the Bethe ansatz out of equilibrium to calculate the 
transport current in the steady state of a quantum dot 
Q . In this work IRLM, particulary, served as a simplest 
test and toy model. The Mehta and Andrei showed that 
the transport current in the interacting resonance level 
model is changed non monotonically with the strength of 
Coulomb coupling. The analysis of this result in refer- 
ence Q was provided by time-ordered scattering formal- 
ism in the next-leading logarithmic order and, indeed, 
it was found a correction to the current in [3J|. Close 
to the work Q a new method to calculate the transport 
current at the first order in Coulomb interaction was sug- 
gested by Doyon [B} . This method takes into account the 
jump of the wave function at the impurity site (impu- 
rity condition) and exploits the Hershfield's approach [f| 
to nonequilibrium dynamic. Our formula for the trans- 
port current slightly deviates from his result and we also 
have obtained the current more directly with the help of 
Keldysh technique. The bosonization of IRSM hamilto- 
nian made in Q was used to calculate linear conductance 
in the region for Coulomb coupling strength where the 
tunnelling remains relevant. The bosonized form of the 
hamiltonian permits a clear identification of this region 
and shows that for sufficiently large Coulomb coupling 
the hybridization energy flows to zero. 

The principal purpose of this work is the calculation 
of the current-current correlation function, in particular, 
the zero frequency shot noise power-the value This value 
is very important experimentally measured characteris- 
tic and technically it is much more difficult to obtain. 
Shot noise has not be considered in the cited articles and 
this value, probably, causes a special difficulties for Bethe 
ansatz approach in nonequilibrium. Here we obtain the 
shot noise power for a general position of bear impurity 



level €d relative to the Fermi energy which is taken at 
zero. 

The hamiltonian: The standard form of the hamilto- 
nian of IRLM has a form 

H = ff +i(X^(°) dt + h.c) + e d dU (1) 

i 

i 

The first term corresponds to non-interacting electrons 
in the two leads 



Ho = ^£j(fc)V'| fc , 0ifc 



i k 

where ipik, £i(k) are the electron field operator and the 
electron energy of a lead i, respectively. Index i = L. R 
indicates left (right) lead. The second term describes the 
tunnelling processes. The last one presents an important 
interacting term with Coulomb capacitive coupling U . 
For non-equilibrium transport properties we use Keldysh 
technique and go to the action generalized for Keldysh 
space. We develop a perturbation theory in U, therefore, 
separate the action A into two parts A = Aq + Ajj where 



A v = -U J *[^^(0) t ^(0)d t d] p ( rf 



(2) 



here o~ z stands for Pauli matrix, the summation over re- 
peated Keldysh indices p = 1,2 is performed. Also [..] p 
means that to every operator in the brackets we have to 
assign Keldysh index p. The remaining part of the action 
is easily diagonalizcd by a shift of electron operators in 
the leads 



tp 



4 P + t(d^ z g ik Y , = # fc + t{g lk o z df (3) 



Here gik are the matrix Keldysh Green functions (GF) of 
the electrons with energy e(k) in the lead i. Actually we 
need GF at the impurity site x = 0: <?i(0) = 27rA(0).gi. 
Here N(0)is the density of states in the wire and gt — 
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/ dekgik- The Keldysh rotation W — ( J ^ 

leads to a especially simple form of these functions: 
g' n = Wa^W- 1 



9* 9« 
9* 



(4) 



where g^ A = =r/2; gf, R {E) = -i taah[{E ± eV/2)/2T] , 
T is the temperature and V -the applied voltage. Below 
we consider the case of T = 0. With the help of ([3]) the 
action Aq then becomes 



Ao = J [ dtdt'^lg^ik + ^G^d} (5) 

ik 

here modified Keldysh GF of impurity level has appeared. 
After Keldysh rotation (QJ this GF transforms to G' 



G'- 1 =G' - 1 -T(g' L +g' R ) 



(6) 



where T = 2irN(0)t 2 stands for hybridization energy and 
Go is the GF of a bear noninteracting level. 

Shot noise: The current noise power is given by 
current-current correlation function Stt 1 — l/2(< 
I(t)I(h) > + < I(h)I(t) >)- < I > 2 which in the 
Keldysh time ordering form can be rewritten as 



St 



l/2(< fl 2 (t)l\t') > + < fl 1 ^)! 2 ^) >) (7) 



Here the superscripts 1,2 are related to the branches of 
the time contour. Current operator is obtained by com- 
mutation of the density operator of the left or right lead 
with the hamiltonian ([2]). In a symmetric presentation 
the current operator acquires a form 



I = 



ite 



<^{Q)U-dH{0)>-^ = ^ L ~^ R (8) 



Next we insert © into formula for noise ([7]) and obtain 
now a working expression to which the Keldysh technique 
is applied. In this formula the average is performing with 
the total action Aq + Ajj. Considering the shift ((3]), we, 
at first, calculate the averaging of equation 1(7} with the 
action Aq. As a result we get transport current Iq = f^i- 
and shot noise So for noninteracting resonance level 



So 
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Me^ + r 2 )- 1 ] 



(9) 
(10) 



where here and below j± = arctan ± arctan -j^- and 
e± = td ± eV/2. In the linear regime (small voltages: 
V/(e 2 + r 2 ) << 1) we immediately recognize standard 
expressions I = (e 2 /h)f, S = {e 3 /h)f{l - T), where 
T = r 2 /(r 2 + e^) is the transmission coefficient. In 
the last formulas we have returned to the standard units 
(here 2ir — > h). 



To consider the Coulomb interacting contribution we 
expand the average in the noise |7J) and current formulaes 
to the first order in Ajj. It is interesting to notice that 
in the Keldysh approach the only divergencies for the 
shot noise and the transport current are of logarithmic 
type which can be accounted by using the renormaliza- 
tion group. Again performing the shift of leads operators 
((3|) we result in a sum of time ordering products for ip^ 
and d? operators which have a simple diagrammatic pre- 
sentation. After rather long, though, direct calculations 
we arrive at the shot noise power S = Sq + Sjj where 



Su 
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or 



uT 

47T 

r 2 [(e 2 , + r 2 )- 1 -(e 2 _ + r^ 
r 2 [(e 2 , + r 2 )- 1 + ( e 2 _ + r 2 )- 1 

ln[ [(e 2 h + r 2 )(e 2 _ + r 2 )]i/4 ] 



- is)} 

(11) 



Similarly, the capacitor coupling leads to a contribution 
to the transport current 



Iu = 2^ {{j -- Jl)L+ 8 lU ^ 



Y 2 1 

p2 + 2 J2J + ^ 



In equations (jlip and P^|) D is the large cut-off of the 
order of band width, and we also introduced the notation 
u = 2irUN(0). It is straightforward to prove by direct 
check that the shot noise power, the same as the current 
H, satisfies the same Callan-Simanzik equation with 0- 
function (3 = —uT/ir 



D dD + ^0T =° 



(13) 



Renormalization group equation for hybridization energy 
follows immediately /x = — where /i is a physical 
cut-off. This equation can be solved in the scaling limit 
by introducing a parameter Tk- This parameter, like 
Kondo temperature, can be defined at the scale /i = D 



and the bear value of T: T, 



1-\-u/tt 



K 



TD U ^. The solution 
gives hybridization coupling as function of physical cut- 
off (U 



r(/x) = 



1 + 7T 

K 



(m) 



-u/n 



(14) 



There is one important remark. Let us use definition of 
Tk in terms of high energy cut-off D and bear tunnelling 
width (as given above). For repulsing interaction (U > 
0) we get an increase of T(/i) and, correspondingly, the 
increase of the average current. This happened because 
D is the largest scale D > /j,. However, in the universal 
region we can properly fix the scale Tk and obtain the 
decrease of the current as U grows. Even at U—0 the shot 
noise is a nonlinear function of hopping probability. In 
this case the Fano factor F = S/I reaches the maximum 
value at T = 1/2, and F = at resonance position of 
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FIG. 1: Fano factor as function of applied voltage for dif- 
ferent values of u. The numerical values of chosen parame- 
ters: u/Tk = 3,€d/Tk — —0.5. The Coulomb coupling u/tv is 
changed from (lowest curve) up to u/n=0.1 



the level (e c i = 0). In Figl. we plot the Fano factor 
as a function of voltage for different values of Coulomb 
coupling U . The lover curve correspond to U = 0. As 
to the physical cut-off parameter /x, it is smaller then D, 
actually ji « D. In the scaling renormalised regime 
moderate change of /i has, indeed, little impact on the 
noise value. 

Strong Coulomb interaction: The repulsive case for 
large U was investigated by renormalization group ap- 
proach [4J with a conclusion that that the orthogonality 
catastrophe reduces the hybridization energy consider- 
ably. Here we use the bosonization of the hamilto- 
nian ([2]) as a method, which is especially convenient for 
application to the present problem. After unfolding and 
linearizing near the Fermi surface, one ends up with the 
hamiltonian Hq for two chiral right movings electrons 
R n (n=l,2). Next, applying the unitary transformation 
Y = exp[ieVt/2N v ] {N v = J dxR^Rj) [| we move 
all dependencies on the applied voltage to the tunnelling 
part of Hamiltonian. Recalling the form of the current 
operator it is naturally to include the current as a source 
part a(t)I to an effective tunnelling hamiltonian 



H T = t[(l + ia)i?i(0)e- ley */ 2 + 

{l-ia)R 2 {0)e leVt f 2 ]S + h.c 



(15) 



We write the hamiltonian in terms of even and odd 
combinations R± = (Ri ± R-i)l\[2. Let us bosonize 
R± = le'^i. This yields 

H = Y^H {ip n ) + ^={d xV ^) + d x ^ + {Q))Sd + 



n=± 
t 



[i 1+ c{t)e lV ^ tp + (Q) +ia(t)ri- 



iy/ 4-inp- 



]rft + h.c 



Here H$(ip n ) = \ J dx(d x ip n ) 2 . The time dependent 
functions c(t),a(i) have a form c(t) — cos(eVt/2) + 



a{t)sin(eVt/2) and a(t) = a(t) cos(eVt/2) - sin(eVt/2). 
If we pass to the action on Keldysh time contour then 
the source function a is transformed into a diagonal 
matrix with the entries a\ and ct2- Next, to cancel 
the direct Coulomb interaction we perform the unitary 
transformation K = e^ d d ( ip - +v +^ \ At this stage, 
we use the new canonical boson variables [7] $± = 
[(V35F- - ^<P T ]/P with p = |( U -7r) 2 + 27r and 

transform the hamiltonian into the desired form conve- 
nient for constructing the Keldysh effective action. 



H = Vff o (0n) + 4=[c(t)»? + e 



The effective action with this hamiltonian can be used 
also in the limit of small U. For this we expand the di- 
mension parameter (3 to first order in U and return to 
fcrmion representation. Here we notice that the whole 
perturbation now is related to the tunnelling hamilto- 
nian. The renormalizcd hopping rate becomes irrelevant 
perturbation if 1 > 87r or \U— tt\ > \^3tt. In this case the 
nonlinear dependence of shot noise power on transmission 
rate is not important and, thus, Fano factor F — ► 1. 

Details of calculations: Our principal equation (III|) as 
well as the equation for the transport current follow as 
a result of summations of a series of Keldysh diagrams 
representing contributions of the first order in U. We get 
general expression for Sjj by inserting the action Ajj and 
the current operators given by eq.© under the average 
in 0. 



-dHmii^d-d^mi, > 



(16) 



where the summations over repeated indices is implied. 
Recalling the shift ([3]) we can express (|I6[) only in terms 
of level d and tilde leads operators. All operators now 
are written in the interaction representation relative the 
quadratic action Aq. After performing transformation 
([3]) shot noise power Su is separated into a three blocks 
having, correspondingly, four, two and none of tilde elec- 
tron operators. As a representative example we consider 
here the simplest black with two pairs of tilde operators 
(like ipi (0)^1(0)). Due to Wick theorem with quadratic 
action Aq the contractions of electron's and impurity d 
operators are separated and we end up with four dia- 
grams presented on Fig. 2. Explicitly, only two diagrams 
(b) give nonzero contribution. We can write this contri- 

U 



bution as S v — ^j-Qu, where 



/ ^{e)gf{e)o-T j ^\e') 9 ^(e') 
(1^2) (17) 
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a) 




b) 

2 Z- V %J 2^ "V j1 

FIG. 2: Keldysh diagrams which contribute to Sp. Solid lines 
represent the leads electron GF <?™'' 3 , the dash ones stand for 
dot GF G(e). The numbers 1,2, s are the Keldysh time in- 
dices. The topologically similar diagrams give, nevertheless, 
different contributions due to rearrangement of Keldysh in- 
dices 



power in the limit of zero temperature in interacting res- 
onance level model. We have shown that in spite of ex- 
isting scepticism against the application of Keldysh the- 
ory to nonequilibrium stationary interacting problems, at 
least, for the weak Coulomb coupling in IRLM this theory 
is justified. In the strong interacting limit we construct 
effective action based on bosonization of the Hamilto- 
nian. We was shown that the orthogonality catastrophe 
causes the Fano factor to approach one. It is important 
that the Coulomb coupling is not renormalized, (U does 
not flow), therefore, renormalized perturbation theory is 
correct for all value e<j, V, Tk -C-D- 

This work was inspired by the lectures given by N. 
Andrei on Bethe ansatz in nonequilibrium. I would like 
to thank B. Horovitz and Y. Meir, for discussion. 



Green's functions G and gi in Eq. fTTj) are related 
to Keldysh GFs g',G' (see (gj|, ©). In particular, 
G R A {e) = (e - e ± zT)" 1 and G K = G R \gf + g%)G A . 
In the limit of zero temperature he energy integrations 
can be easily performed with the result 

Qu = -^Lj_ (18) 

The other blocks of diagrams include more functions gi 
and G, however, the energy integrations are simple and 
can be easily done. 

Conclusion: We used standard Keldysh technique to 
calculate for the first time zero frequency shot noise 
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